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For low-dimensional metallic structures, such as nanotubes, the exchange coupling between lo¬ 
calized magnetic dopants is predicted to decay slowly with separation. The long-range character 
of this interaction plays a significant role in determining the magnetic order of the system. It has 
previously been shown that the interaction range depends on the conformation of the magnetic 
dopants in both graphene and nanotubes. Here we examine the RKKY interaction in carbon nan¬ 
otubes in the presence of uniaxial strain for a range of different impurity configurations. We show 
that strain is capable of amplifying or attenuating the RKKY interaction, significantly increasing 
certain interaction ranges, and acting as a switch: effectively turning on or off the interaction. 

We argue that uniaxial strain can be employed to significantly manipulate magnetic interactions 
in carbon nanotubes, allowing an interplay between mechanical and magnetic properties in future 
spintronic devices. We also examine the dimensional relationship between graphene and nanotubes 
with regards to the decay rate of the RKKY interaction. 


I. INTRODUCTION 

Recently, there has been much investment in the 
field of spintronics motivated by the tremendous po¬ 
tential for technological applications. Low-dimensional 
structures such as graphene^, nanowires^, nanotubes^, 
nanoribbons^, silicene® (and many more) are expected 
to lead to useful spintronic applications, possibly lead¬ 
ing to the production of extremely efficient magnetic 
sensors, high-capacity memory storage, and non-volatile 
computer memories®’^. Important to spintronics is the 
mechanism of interaction between embedded impurities 
known as the indirect exchange interaction. 

This interaction is one of many such interactions me¬ 
diated by the conduction electrons of the host mate¬ 
rial and is realized as the energetically favourable con¬ 
figurations of localized moments, driven by the energy 
difference between them. Usually calculated within 
the Ruderman-Kittel-Kasuya-Yosida (RKKY) approxi- 
mation®^^^, this interaction has been extensively studied 
in graphene^^^^^ and carbon nanotubes^^^^®, where the 
focus has been on the sign, magnitude, and rate of de¬ 
cay of the interaction^"^’®^^^^. This type of analysis has 
been carried out for the 4 main types of impurities: sub¬ 
stitutional, top-adsorbed, bridge-adsorbed, and center- 
adsorbed (Fig. 1). These impurities are differentiated by 
their conformation with the lattice, which is found to 
have a strong effect on the behaviour of the RKKY in- 
teraction^®’^^’^^’^®’®®’"^^. In carbon nanotubes (CNTs) the 
RKKY interaction is thought to decay as D~^, where D 
is the separation between impurities. This is the case for 
substitutional, top- and bridge-adsorbed impurities, im¬ 
plying that the magnetic moments of adatoms are able 
to feel their mutual presence even when they are very 
far apart. For center-adsorbed impurities (also known 
as plaquette impurities) the interaction is predicted to 
have a D~^ decay rate®®. This decay rate is significantly 


faster than for other impurity types, and is not as useful 
for applications. This is unfortunate as it is the preferred 
configuration for many common Transition Metals, which 
are likely magnetic dopants^®^"^®. 

Controlling the interaction between magnetic objects 
in carbon nanotubes (CNTs) may provide an entry into 
the technologically promising area of spintronics. In or¬ 
der to expand the applicability of those systems we need 
to understand how the interactions may be modified to 
suit our needs. Much progress has been made explain¬ 
ing the intrinsic properties of nanotubes and nearly two 
decades after their discovery CNTs are still the subject of 
intensive scientific research due to their intriguing phys¬ 
ical properties'^^^^®. One avenue that has been explored 
to alter the intrinsic properties of nanotubes is the in¬ 
troduction of strain. The effects of strain on the elas¬ 
tic, structural, and electronic properties of carbon nan¬ 
otubes - including band gap and electron conductance - 
have been extensively studied®®^®®. Many uses have been 
suggested for strained nanotubes including strain sensors 
and Field-Effect Transistor®^’®®. It has been shown that 
uniaxial strain can open a bandgap in metallic nanotubes 
and alter the bandgap in semi-conducting nanotubes®®’®®. 
One avenue that has not been explored is how strain may 
affect their magnetic properties. 

It has been shown that the properties of RKKY in¬ 
teraction in both graphene^^’®®’^® and nanotubes®®’®® are 
strongly influenced by the conformation of the magnetic 
impurities with the lattice. In graphene it has been 
shown that uniaxial strain has the ability to amplify or 
attenuate that interaction, as well massively increase the 
range or suppress the interaction in some cases®®’^®’®^’®®. 
In graphene nanoribbons strain has been shown to tune 
the exchange splitting of non-vanishing moments induced 
by vacancies in the lattice®®. In this paper we determine 
the role that strain plays in controlling the RKKY inter¬ 
action in achiral carbon nanotubes, and explore the pos- 
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FIG. 1. A schematic of an armchair nanotube with impurities, 
represented by the red circle, whose interactions with the nan¬ 
otube lattice are indicated by the solid green lines. Four differ¬ 
ent types of impurities are represented: a) substitutional im¬ 
purities, b) top-adsorbed impurities, c) bridge-adsorbed im¬ 
purities, d) center-adsorbed impurities, e) A schematic of the 
unit cell, highlighted in blue, and the inter-atomic hopping 
integrals t 2 and t\. 


sibility of using strain to modify the interaction between 
such impurities. We examine the four main types of 
impurity conformations: substitutional, where a carbon 
atom is replaced by an impurity; top-adsorbed, where 
an impurity sits above, and connects to, a single car¬ 
bon atom; bridge-adsorbed, where an impurity sits above 
the bond between two neighbouring carbon atoms; and 
center-adsorbed, where an impurity sits at the center of a 
carbon hexagon and attaches equally to all six surround¬ 
ing carbons. This work is carried out for both armchair- 
edged nanotubes (ACNT) and achiral zigzag-edged nan¬ 
otubes (ZZNT) as they have substantially different re¬ 
actions to strain. Figure 1 schematically shows the four 


conformation types in an achiral ACNT, where the impu¬ 
rities are separated along the axial direction of the tube. 
We provide simple expressions for the interactions be¬ 
tween two impurities on a strained nanotube. The sim¬ 
plicity of these expressions allow us to understand the 
modified coupling in terms of the inter-carbon hopping 
integrals, and analytically understand how strain may 
change the decay rate drastically. We also explore the 
role played by nanotube circumference in determining the 
decay rate of the RKKY interaction between impurities 
in order to validate our model. 


II. METHODS 

The energy difference between the ferromagnetic (FM) 
and antiferromagnetic (AFM) alignments of two mo¬ 
ments embedded in a conducting host is described as the 
indirect exchange coupling between the two moments. 
The Lloyd formula method allows us to calculate the to¬ 
tal energy difference. Jab, between two magnetic impu¬ 
rities labelled A and B. The Lloyd formula is given by 

Jab = -^JdE f{E) In (l + 4^^ gab{E)gUE)) (1) 

where f{E) is the Fermi function, I4x is a spin-dependent 
onsite potential that accounts for the exchange splitting 
in the magnetic orbitals, and GAsiE) is the real-space, 
single-electron Green Function (GF) describing the prop¬ 
agation of electrons with spin a =t or J,. According 
to our definition J < 0 (J >0) corresponds to a fer¬ 
romagnetic (anti-ferromagnetic) alignment of magnetic 
moments. Our nanotube is modelled using the Nearest 
Neighbour Tight Binding Approximation with an inter¬ 
atomic hopping integral to = —2.7el4, which provides a 
good approximation of the electronic structure of GNTs 
for all but the smallest circumference tubes. Unless oth¬ 
erwise specified we will use to as our energy unit. 

We assume that each substitutional impurity causes a 
change in onsite energy, and that each adsorbed impu¬ 
rity orbital has a finite hopping r to adjacent sites 
on the graphene lattice: top-adsorbed {N = 1), bridge- 
adsorbed {N = 2), and center-adsorbed {N = 6). We 
will only consider a single magnetic orbital at each im¬ 
purity site, separated along the longitudinal direction of 
the nanotube, however, it is straightforward to generalize 
this approach to deal with multiple orbitals, or separa¬ 
tions with an axial component. We will not deal with the 
exact parametrizations for specific impurity types here. 
These parametrizations can be found in numerous ab ini¬ 
tio studies^^’^^’^® of single impurities in nanotubes with 
different configurations, though it is worth noting that 
the results presented here are not strongly dependent on 
the impurity type used. 

Strain is introduced into the nanotube system, in a 
similar manner as for a graphene system^®^^’®"^’®®, by 
splitting the hopping integral between carbon atoms, to) 
into t 2 and ti. In the convention used here is the 
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iiitra-unit-cell hopping and ti is the inter-unit-cell hop¬ 
ping (Fig. le). The strain-dependent, real-space GF be¬ 
tween two sites on the graphene lattice can be written as 
a double integral over the Brillouin zone^®’^^ in the form 


1 

gab - ^ 


7r/2 TT 
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dkA 


N^{E,k) 

E^-\f{kW ’ 


( 2 ) 


where 


f{k) =t2 + 2tiCOs{kz)e'’^'^. (3) 

Here contains information about the sites of each im¬ 
purity, D is a separation vector, and /(k) is a sum of 
Bloch phase terms over nearest neighbours. 

The Green’s function (GF) for nanotubes is derived in 
a similar manner to that of the GF for pristine graphene, 
except that the periodicity in the nanotube causes one of 
the components of one of the k-vectors to be quantized. 
The integral then becomes a sum over all the unique k- 
points, which are obtained from the quantization condi¬ 
tion which describes the circumference of the nanotube, 
ric- Following the procedure of Ref. 32 the GFs for both 
ACNTs and ZZNTs may now be written in a general form 

(4) 
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which will allow us to solve the Lloyd formula analyti¬ 
cally. 

To calculate the coupling between two impurities, A 
and H, we now use the RKKY approximation to write 
Eq. 1 as 

J^lvaj dEf{E)Vlg 

AB- (5) 

Using the general form of our GFs we can simplify our 
coupling to 


III. RKKY INTERACTION AS A FUNCTION 
OF NANOTUBE CIRCUMFERENCE 

It has previously been reported that in nanotubes the 
RKKY interaction between substitutional, top-adsorbed, 
and bridge-adsorbed impurities decays as with a = 
1, and a = 5 for center-adsorbed impurities, where D is 
the separation between the impurities^^. However these 
have been calculated without much regard for the effect 
that the circumference of the nanotube has on the in¬ 
teraction. In bulk graphene these rates are found to 
be a = 3 for substitutional, top-adsorbed, and bridge- 
adsorbed, and a = 7 for center-adsorbed^^. In this sec¬ 
tion we examine the role of the nanotube circumference, 
Uc (which can be easily related to the more usual measure 
of diameter), for both armchair nanotubes and metallic 
zigzag nanotubes, to understand how the RKKY inter¬ 
action transitions from the behaviour predicted in nan¬ 
otubes to that predicted in bulk graphene. Here ric is 
a dimensionless quantity which describes the number of 
vectors between unit cells, in the armchair (zigzag) direc¬ 
tion for AGNTs (ZZNTs), that are required to traverse 
the circumference of the nanotube. 



J ^ ^ Im/ dEf{E)VlBiE,j, (6) 

3,k 

where 6(A, j, fc) = A{E, j)A{E, k). 

We have shown previously that the magnetic coupling 
can be easily extracted when the coupling is expressed 
in such a form, by reducing the integration to a sum 
over Matsubara frequencies and expanding the functions 
B{E) and Q{E) around the Fermi energy in the low tem¬ 
perature limit^®’^°. Note that B{E) and Q{E) are also 
strain-dependent, but the variable, e, has been omitted 
for neatness. The coupling between impurities can then 
be expressed as 


j - X! 

j,k 


^1iQ^°'>D 




S(0 


_(2iQ(i))Ui Dd+^)_ 


(7) 

where the dependence of Q and B on A, j, k have been 
omitted for simplicity. 


FIG. 2. Decay (o) versus circumference (ric) for ACNT (solid- 
black) and ZZNT (dashed-blue) for a) substitutional impuri¬ 
ties, same-sublattice, b) substitutional impurities, opposite- 
sublattice, c) bridge-adsorbed, d) center-adsorbed. 


Figure 2 shows the change in decay, a, between impu¬ 
rities as a function of the circumference for both ACNT 
and ZZNT at Ep = 0. For both same- and opposite- lat¬ 
tice substitutional impurities, and bridge-adsorbed im¬ 
purities we see that the greatest predictor of behaviour 
is whether the nanotube is an ACNT or a ZZNT. On AC¬ 
NTs the interaction very quickly decays away from the 
nanotube case of a = 1 to the graphene case of a = 3. 
This is seen in 2 a), b), c) as the black line, which re¬ 
sembles graphene by about ric = 200. On ZZNTs the 
interaction has a plateau of NT behaviour before mov¬ 
ing gradually to the behaviour seen in graphene, this is 
shown as the blue line, and does not resemble graphene 
until about Uc = 1000. 
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Figure 2 d) shows the same effect for center-adsorbed 
impurities. By increasing the circumference of the NT 
we show that both ACNT and ZZNT decay rates go from 
to Again, impurities on ACNT decay much 

more rapidly towards the bulk case and on ZZNT exhibit 
a plateau and a slow decay towards bulk. In ACNTs the 
axial separation is along the zigzag direction, which is 
known to feature a repeating triplet pattern of interac¬ 
tion strength: strong, strong, weak^^. The circumfer¬ 
ence was increased in steps of three to produce a smooth 
curve. This result shows a geometry dependence in the 
transition from nanotube-like to graphene-like behaviour 
as the diameter of the CNT increases, and may be rele¬ 
vant for studies of the largest feasible CNTs since there 
is a notable difference between the two geometries even 
at quite small tic- 


IV. AXIAL STRAIN EFFECTS ON RKKY 
INTERACTION IN NANOTUBES 


£ 



E 


FIG. 3. The density of states, p, at E = 0 for a ZZNT with 
ric = 21 under uniaxial strain, e. Subfigures show the density 
of states against energy for three strain values: e = 0 (red), 
e = 0.05 (blue), and e « 0.11 (green). 


Unstrained achiral nanotubes may be metallic or semi¬ 
conducting. ACNTs are known to be metallic for all 
circumferences, ric, whereas ZZNTs are only metallic if 
Tie = 3k, where k is an integer, and are otherwise semi¬ 
conducting. These properties are due to the intersection 
of the discretized momentum values and the Fermi sur¬ 
face. The introduction of strain, e, into the nanotube 
system alters the coincidence of the momentum values 
and the Fermi surface, this in turn changes the electronic 
properties of nanotubes. ACNTs remain metallic under 
strain, however strain can change a conducting ZZNT 
into a insulating one, and conversely an insulating one 
into a conducting one^®. Figure 3 shows the density of 
states (DOS) of a ZZNT with Uc = 21 as strain is ap¬ 
plied. The initially metallic NT immediately becomes 
semi-conducting with slight strain, and as strain is fur¬ 
ther applied the two peaks in density of states come to¬ 
gether and eventually merge. This is seen in all circum¬ 
ferences of semi-conducting ZZNTs, where the precise 
value of strain required to shift the density of states away 
from zero around Ep = 0 decreases as the circumference 
increases. 


A. Armchair-edged Nanotubes 

Although ACNTs remain metallic with strain, this 
does not mean that the RKKY interaction is unchanged. 
Figure 4 shows the effect of strain on the RKKY interac¬ 
tion for substitutional and bridge-adsorbed impurities for 
several separations D, on an ACNT with ric = 21. Top- 
adsorbed impurities have been omitted due to the simi¬ 
larity of behaviour with substitutional impurities. In all 
three cases oscillations are seen, but only in c) are these 
oscillations around zero signifying a change from ferro¬ 
magnetic to anti-ferromagnetic behaviour. For both a) 
same-sublattice and b) opposite-sublattice cases, strain 


is seen to either amplify or attenuate the interaction de¬ 
pending on the precise value of strain. This allows for the 
magnitude of these interactions to be precisely controlled. 
Interestingly there exist strain values (well within the ex¬ 
perimental limits) that completely shut off the opposite- 
sublattice interaction, whilst simultaneously maximizing 
the same-sublattice interaction. In addition to this strain 
seems to generally amplify the interaction, while leaving 
the decay rate unchanged. This is understood from Eq. 7. 
The overall decay rate can be found from the first non¬ 
vanishing S(e), which, for same-sublattice substitutional 
impurities, is 




I 

iric'^itl - itj)' 


( 8 ) 


This factor corresponds to a coupling decay rate of a = 1, 
identical to that of the unstrained case. The equivalent 
opposite-sublattice term is similar. Small values of strain 
will not cause this term to vanish, and so no change in 
the decay rate is expected. 

For bridge-adsorbed impurities (Fig. 4 c) strain is ca¬ 
pable of switching the sign of an interaction in a very 
controlled fashion. Since bridge-adsorbed impurities dis¬ 
play a coupling with alternating sign as a function of 
separation strain cannot create wholly FM or AFM in¬ 
teractions between multiple impurities. However, for a 
pair of impurities it is capable of switching the sign of 
the interaction. The coupling between bridge-adsorbed 
impurities is weakened with strain, however the decay 
rate is not affected. The first non-vanishing B{e) is 


B{e)^°^ = 


4 

ric'^itl - 4tf) ’ 


(9) 


again corresponding to a decay rate of D ^ for small 
strains. 
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FIG. 4. Plot of the coupling (J) as a function of strain (e) 
for ACNTs at D = 10 (black) and D = 20 (dashed-red) for a) 
same-sublattice substitutional, b) opposite-sublattice substi¬ 
tutional, and c) bridge-adsorbed impurities, where the dotted- 
black line indicates the transition from FM to AFM. 
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FIG. 5. The interaction between two centre-adsorbed impu¬ 
rities in an AGNT in a) the absence of strain (e = 0), and 
b) the presence of strain (e = 0.1). The inset is a normalized 
log-log plot of both interactions. The coloured lines show the 
decay rates: (dash-dotted-blue), D~^ (dashed-red). 


For center-adsorbed impurities we see that the intro¬ 
duction of strain has a large impact on the decay rate and 
hence the magnitude of the interaction. Figure 5 shows 
the behaviour of the interaction between center-adsorbed 
impurities with and without strain. Strain changes the 
interaction from being primarily antiferromagnetic to be¬ 
ing split evenly between ferromagnetic and antiferromag¬ 
netic. The interaction between these impurities is known 
to decay rather quickly as D~^ (Fig. 5 a). This can be 
understood from the terms, which are zero for 

Z = 0,1, 2,3. The first non-zero term is found at 
(it is too long to reproduce here) and thus indicates the 
known decay rate of D~^. The small values here owe to 
the necessity of examining the coupling at large ranges 
to ascertain a decay rate. 

However, the introduction of strain breaks a key sym¬ 
metry of the system massively increasing the range of the 
interaction. The interaction in the presence of strain de¬ 
cays slowly as D~^ (Fig. 5 b). The inset of figure 5 shows 
a log-log plot of the coupling between the impurities for 
the unstrained and strained cases. Examining the 
we find an explanation for the increased range. 


.( 0 ) ^ (t2-^l)V(£) 
AncH\{tl - 4£l) 


( 10 ) 


where is a strain-dependent phase term. This simple 
expression allows us to understand the role that strain 
plays on the decay rate of centre-adsorbed impurities 
Here the numerator displays a. t 2 — ti term, since the 
unstrained system is defined by t 2 = ti it is exactly zero 
only in unstrained graphene. Such a non-zero B{e)^^'^ 
term corresponds to a decay rate of D~^. This symme¬ 
try breaking effect is analogous to that seen in graphene 
where the decay rate is massively reduced from D~^ to 
D~^ by the introduction of strain^°. This effect allows 
the strength of the interaction to be massively increased 


with a relatively small amount of strain. 


B. Zigzag-edged Nanotubes 


We now move our attention to ZZNTs. Since the in¬ 
troduction of strain changes a conducting ZZNT to a 
semi-conducting ZZNT the effect of strain in ZZNTs is 
markedly different. Figure 6 shows the coupling between 
substitutional, bridge-adsorbed, and center-adsorbed im¬ 
purities in ZZNT with = 21 as a function of strain. 
Since the introduction of strain changes a ZZNT from 
conducting to seini-conducting the RKKY interaction 
in an initially conducting ZZNT goes to zero (Fig. 6 
a). Strain also has the opposite effect on initially semi¬ 
conducting ZZNTs, taking the RKKY away from zero 
for precise values of strain. Similar behaviour is seen for 
bridge-adsorbed impurities (Fig. 6 b), with the key differ¬ 
ence being the sign change that occurs around this precise 
strain value. What this means is that for two bridge- 
adsorbed impurities separated by some distance D there 
exists a small range of strain that can quickly turn the 
RKKY interaction from strongly AFM to strongly FM. 
This is in contrast to ACNTs where a sign change can be 
achieved but only gradually. The precise strain value will 
depend on the circumference of the nanotube ric- In ad¬ 
dition, the sign change occurs at all separations, meaning 
that it is possible for strain to quickly change a system 
of N impurities from all having the same alignment to 
having alternating alignments leading to no overall mag¬ 
netic moment. The center-absorbed impurity (Fig. 6 c) 
shows similar behaviour, however since a key symmetry 
that leads to the D~^ decay is broken by strain the mag¬ 
nitude is increased massively, which gives the appearance 
of initially being zero. This could prove extremely useful 
for spintronic applications where it allows minor changes 
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FIG. 6. The coupling between a) substitutional, b) bridge- 
adsorbed, and c) centre-adsorbed impurities is plotted as a 
function of strain in a ZZNT with ric = 21. Each coupling 
is plotted for three values of separation D = 10 (solid-black), 
Z) = 20 (dashed-red), D = 30 (dash-dotted-blue). The grey 
bands indicate the regions of sudden FM/AFM switch. 

in strain to result in a large change in the magnetic or¬ 
dering of impurities. This change could be implemented 
in a reversible and controllable manner. 


V. CONCLUSIONS 

In summary, we have presented simple expressions for 
the RKKY interaction between magnetic impurities in 
strained nanotubes. We first validated our model by 
examining the role that the circumference of the nan¬ 
otube plays in the decay rate as the nanotubes become 
more graphene-like, and found that ZZNT maintain their 
nanotube-like behaviour for larger circumferences than 
ACNT. We have also shown that strain can change the 
magnitude, sign, and decay of these interactions in AC- 
NTs. The amplification was particularly pronounced for 


centre-adsorbed impurities where we showed that the 
symmetry breaking of the hexagonal lattice by uniax¬ 
ial strain leads to a significantly slower decay rate be¬ 
tween centre-adsorbed impurities: D~^ to D~^ in arm¬ 
chair nanotubes. Zigzag nanotubes, meanwhile, display 
a wide range of amplification and switching effects with 
minor variations of the applied strain. These features are 
related to the transition from metallic to semiconducting 
behaviour. Experiments to date searching for magnetism 
in disordered graphene seem to suggest paramagnetic, 
non-interacting moments®®. Signatures of indirect ex¬ 
change interactions between such moments in graphene 
are very difficult to detect due to their short-ranged na¬ 
ture, particularly if they adopt certain adsorption con¬ 
figurations. Amplification of these couplings using strain 
may provide a path to their detection in future experi¬ 
ments. The abrupt changes in interaction strength and 
behaviour, demonstrated here for ZZNTs, suggest that 
strained nanotubes can act as switches for magnetoresis¬ 
tive behaviour, which may find use in future spintronic 
devices. Thus strain presents itself as a way to control in¬ 
teractions between magnetic dopants in nanotubes, and 
the overall magnetic ordering of the material, in a re¬ 
versible and controllable manner - an indispensable trait 
for future spintronic devices. 
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